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I.  INTRODUCTION 


The  resistance  to  inelastic  deformation  of  amorphous  metallic  and  polymeric 
materials  depends  on  the  state  of  aging  of  their  structure,  i.e.,  the  degree  of  com¬ 
pletion  of  structural  relaxation.  Such  structural  relaxations  have  been  widely 
studied  experimentally  by  many  techniques,  which  include:  direct  probing  using 
X-ray  scattering  [l],  decay  of  excess  enthalpy  by  calorimetry  and  related  means 
[2],  and  indirect  probing  through  the  effect  of  aging  on  mechanical  properties, 
such  as  internal  friction  [3-5],  viscous  flow  [6],  micro  hardness  [4,7],  in  both  amor¬ 
phous  metals,  and  in  polymers  [8].  While  such  structural  relaxations  have  been 
characterized  quite  well  through  changing  structural  parameters,  such  as  radial 
distribution  functions  (RDF),  etc.,  the  actual  kinematics  of  the  local  relaxation 
processes,  which  produce  both  topological  and  chemical  short  range  order,  are 
not  well  understood.  Not  much  is  known  about  how  these  processes  are  related 
to  the  experimentally  observed  kinetics  of  relaxation.  Since  inelastic  deformation 
in  amorphous  media  preferentially  selects  the  portion  of  the  volume  of  excess 
properties,  such  kinematical  features  of  the  relaxation  are  of  primary  interest. 
Because  these  processes  are  complex,  inhomogeneous,  and  very  local,  computer 
simulations  using  the  molecular  dynamics  (MD)  or  the  Monte  Carlo  (MC)  ap¬ 
proaches  should  be  promising  for  their  visualization.  In  the  past,  Takeuchi  and 
coworkers  [9,10]  have  performed  such  simulations  for  CusZri_s  type  glasses  in 
three-dimensional  models  and  have  reported  that  the  principal  driving  force  for 
the  relaxations  is  the  excess  in  atomic  level 
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deviatoric  stresses,  and  that  the  principal  result  is  the  development  of  chemical 
short  range  order.  Kobayashi  and  Takeuchi  [10]  have  also  given  additional  infor¬ 
mation  on  the  changing  distributions  of  atomic  level,  mean  normal  stress,  and 
some  atomic  displacement  fields  that  accompany  the  structural  relaxation  pro¬ 
cess.  In  spite  of  these  useful  observations,  both  the  kinematics  of  local  rearrange¬ 
ments  and  associated  kinetics  remain  inadequately  understood  in  relating  the 
changing  structure  to  its  deformation  resistance  through  local  mechanisms.  We 
have  now  carried  out  a  more  complete  MD  simulation  using  the  same  truncated 
Cu  —  Zr  pair  potentials  used  by  Kobayashi  et  al.  [ll],  but  in  a  two-dimensional 
cell  using  periodic  boundary  conditions  to  better  visualize  the  topological  fea¬ 
tures  of  the  structure.  In  a  preceding  communication  [12],  to  be  referred  to 
here  as  (I),  we  reported  results  of  a  simulation  of  the  melting  and  glass  transi¬ 
tion  process  in  a  two-dimensional  ideal  material  held  together  by  the  same  pair 
potential.  In  this,  and  a  following  communication  [13],  to  be  referred  to  here 
as  (III),  we  report  respectively  the  topological  features  and  the  details  of  the 
kinetics  of  structural  relaxations. 

II.  DETAILS  OF  THE  SIMULATION 

2.1  The  molecular  Dynamics  Method 

The  details  of  the  molecular  dynamics  method  that  was  used  in  the  present 
simulation  were  described  in  more  detail  in  (I).  There,  it  was  shown  that  single 
component  materials  axe  too  unstable  in  the  amorphous  solid,  and 
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crystallize  exceedingly  rapidly,  but  that  two-component  solids  with  appropriate 
size  ratios  such  as  Cu  and  Zr,  are  considerably  more  stable  and  do  not  crystallize 
as  rapidly  when  relaxed  at  low  temperature,  making  them  more  suitable  for 
studies  of  structural  relaxations.  Thus,  the  present  simulation  uses  the  same 
two-dimensional  cell  under  periodic  boundary  conditions  of  a  Cu0.s£r0.s  type 
solid  held  together  with  a  4-8  type  truncated  Lennard-Jones  potential,  as  was 
described  in  (I).  There,  it  was  shown  that  the  fundamental  time  constant  of 
the  atomic  period  of  this  material  was  rj  =  ( mrl/E0)o  s  =  5.40x10“  18sec.,  and 
that  under  an  external  normalized  pressure  of  p*  =  l.O^*)  the  material  melts 
at  a  normalized  temperature  of  T £  =  0.25  and  undergoes  a  glass  transition 
at  a  normalized  temperature  of  0.18.  The  time  steps  between  different  states 
in  the  integration  of  the  equations  of  motion  were  0.02  or  l.lxl0“14sec.  As  in 
(I),  in  the  present  simulation,  the  magnitude  of  the  so-called  border  mass  used 
to  stabilize  the  solutions  and  maintain  the  external  pressure  constant  was  4. 
(See  Appendix  I  of  the  accompanying  communication  (III)  on  the  kinetics  of 
structural  relaxation  [13]  on  the  choice  of  border  mass.) 


(*)For  nomenclature  and  definition  of  normalization  terms,  see  the  Appendix 
of  (I). 
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The  amorphous  medium  and  its  state  need  to  be  described  by  a  set  of 
parameters.  The  ones  that  have  been  found  most  useful  are:  (a)  the  volume  per 
atom  fl,  delineated  in  2-D  models  readily  by  the  Voronoi  polygons  constructed 
around  each  atom  site,  and  the  free  volume,  s/,  which  is  the  excess  of  a  local 
average  in  this  volume  over  and  above  a  certain  reference  volume  for  the  entire 
field  of  polygons,  (b)  the  atomic  level  stress  tensor,  craPi  related  directly  to  the 
distribution  of  material  misfit  in  the  amorphous  material  [14],  (c)  the  atomic 
level  elastic  constants  Ca^j,  related  to  the  state  of  dilatation  and  shear  of  the 
local  site,  (d)  the  atomic  level  enthalpy,  h,  useful  in  describing  the  excess  prop¬ 
erties  of  an  atomic  site,  and  (e)  the  atomic  site  distortion  parameter,  to,  giving 
a  measure  of  the  distortion  of  the  local  atomic  site  compared  to  an  equiaxed  site 
among  uniform  Voronoi  polygons.  We  now  define  these  parameters  below. 

2.2.1  The  Voronoi  polygons  and  free  volume 


In  the  2-D  material,  which  is  the  subject  of  this  simulation,  the  volume  per 
atom  is  represented  by  the  area  delineated  by  a  polygon  constructed  around  an 
atomic  site  by  the  perpendicular  lines  bisecting  the  lines  connecting  atom  cen¬ 
ters.  In  a  perfect  2-D  hexagonal  material,  all  Voronoi  polygons  are  hexagons. 
In  (I),  however,  we  showed  that  when  an  initially  hexagonal  2-D  material 
undergoes  melting,  many  of  the  hexagonal  polygons  are  transformed  into  pen- 
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tagons  and  heptagons,  and  that  these  two  associate  into  structural  dipoles  of  5 
and  7  sided  polygon  pairs,  which  in  turn  associate  into  strings  that  percolate 
through  the  simulation  cell  in  the  melt.  While  these  alterations  in  the  number  of 
edges  of  the  Voronoi  polygons  are  the  most  dramatic  manifestations  of  change, 
other  corollary  changes  also  occur  in  the  number  of  edges  meeting  at  corners  of 
polygons,  i.e.,  instead  of  only  corners  where  3  edges  meet,  other  corners,  where 
4  edges  meet  also,  appear  with  considerable  regularity.  The  meaning  of  these 
special  structural  features  and  how  they  re-distribute  themselves  in  the  melt,  or 
how  they  age  out  at  lower  temperatures,  is  the  main  subject  of  discussion  in  this 
communication. 

Since  significant  increases  in  local  dilatation  results  in  substantial  reductions 
in  cohesive  interaction,  which  has  important  consequences  in  the  deformation  of 
an  amorphous  solid,  the  excess  in  the  local  volume  per  atom  over  and  above  that 
in  a  reference  material,  is  an  important  parameter  of  state.  This  local  structure 
parameter,  which  is  called  the  free  volume,  was  evaluated  in  the  two-dimensional 
model  of  the  present  simulation  as 


1  n 

3/(0  =  -z23<-a°  » 

n  ; 


(1) 


where  a,-  is  the  area  of  a  Voronoi  polygon,  and  the  sum  n  is  carried  out  over 
the  central  atom  and  its  nearest  neighbors  that  enter  into  the  definition  of  its 
Voronoi  polygon  faces.  The  reference  area  a,  was  taken  as  the  average  Voronoi 
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polygon  area  in  the  entire  simulation  ceil. 


2.2.2.  The  atomic  level  stress  tensor 


Egami  and  Vitek  [14]  have  emphasized  the  utility  of  the  atomic  level  stress 
tensor,  defined  much  earlier  by  Born  and  Huang  [15],  to  describe  the  large  local 
variations  of  interatomic  forces  that  result  from  the  distributed  material  misfit 
that  accompanies  the  amorphous  state.  The  atomic  level  stress  is  defined  as: 


*«/»(*)  =  5n~£ 

Mi  jjti 


dr 


r?/ 

1 1 ) 1 X  J 

1 

ra 


(2) 


where  $  is  the  truncated  interatomic  pair  potential  of  the  4-8  Lennard-Jones 
type,  introduced  in  (I),  r<y  is  the  magnitude  of  the  radius  vector  connecting 
atoms  i  and  j,  r"  is  the  magnitude  of  the  component  in  the  local  a  direction  of 
the  vector  connecting  atoms  i  and  j,  and  fl,  is  the  volume  of  the  central  atom 
occupying  site  for  which  the  component  of  the  stress  tensor  is  being  calculated. 
The  sum  is  over  all  near  neighbors  j  surrounding  the  atom  t,  which  are  within 
the  potential  cutoff  radius.  In  the  present  simulation  here,  as  well  as  those  in 
(III)  and  (IV),  what  will  be  of  interest  is  not  only  the  shear  stress  a,„(»)  referred 
to  the  principal  axes  of  the  undisturbed  simulation  cell,  but  also  the  maximum 
shear  stress  r,-  in  the  plane  regardless  of  direction,  and  the  atomic  level  pressure 
Pi,  both  defined  as: 
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Pi  =  " £  (*“(*) +<fw(*)) 


> 


2.2.3  The  atomic  level  elastic  constants 


In  the  amorphous  medium,  where  large  local  variations  in  free  volume  exist, 
forcing  the  material  into  non-linear  modes  of  response,  there  are  important 
variations  in  the  local  elastic  properties  in  the  neighborhoods  of  each  atomic 
site.  These  local  moduli  are  characterized  by  an  elastic  constant  tensor,  which 
is  defined  as  [16]: 


Ca0-7«(O  — 


d29 

diaffdtff 


r  dr 


^7 


dr 2 


189 
r  dr 


> 


(5) 


where  6ai  is  the  Kronecker  delta,  and  the  sum  on  j  is  again  over  ail  nearest 
neighbors,  over  which  the  truncation  permits  interactions.  In  Eqn.  (5),  tap 
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and  e^e  are  tensor  strain  components,  which  are  given  only  to  define  the  elastic 
constants.  They  are  not  explicitly  used  here,  and  will  therefore  not  be  developed. 
In  (IV),  where  imposed  deformations  are  of  interest,  strains  will  be  defined  for 
use  there  [12]. 

2.2.4.  Atomic  level  enthalpy 


The  excess  enthalpy  has  long  been  identified  to  be  an  important  character¬ 
izing  parameter  of  amorphous  materials,  where  it  has  been  considered  by  many 
investigators  as  the  principal  driving  force  in  structural  relaxation.  We  define  it 
as: 


MO  =  \  £  +  pn  •  (6) 


2.2.5.  The  atomic  site  distortion  parameters 

Many  Voronoi  polygons  in  the  melt  and  also  in  the  unrelaxed  glassy  state 
are  considerably  distorted.  These  distortions  are  geometrical  complements  to 
the  free  volume.  Although  both  evoke  atomic  level  internal  stresses  and  changes 
in  elastic  constants,  their  quantitative  description  provides  additional  visual 
evidence  of  the  strong  interactions  in  the  “liquid-like"  material  between  polygons 
which  combine  excess  free  volume  and  distortions.  We  define  the  site  distortion 
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parameter  w(i)  as  the  ratio  of  the  perimeter  length  of  the  actual  polygon  to 
the  circumference  of  a  circle  having  the  same  area  as  the  polygon,  i.e., 


»w  = 


(7) 


where  s,-  is  the  actual  area  of  the  polygon.  Thus,  for  a  circle  w(t)  =  1.0  by 
definition.  This  is  the  minimum  value  for  any  u>(t),  and  corresponds  to  the 
highest  degree  of  symmetry  that  a  polygon  can  have. 

The  site  distortion  parameter  w(i)  increases  with  decreasing  numbers  of 
edges  of  a  polygon  for  regular  polygons.  This  is  a  trivial  increase,  which  par¬ 
tially  dilutes  the  more  relevant  and  interesting  increase  that  comes  form  true 
distortion.  These  trivial  base  level  increases  in  w„  for  regular  polygons  are  listed 
in  Table  I. 

More  importantly,  the  distortion  parameter  increases  with  increasing  depar¬ 
ture  from  regularity  of  a  polygon.  To  separate  the  real  site  distortion  from  the 
background  increases  listed  in  Table  I,  it  is  useful  to  define  a  net  site  distortion 
parameter  tp'(t),  which  gives  the  difference  between  the  actual  measured  valued 
and  the  background  value  defined  as: 

w'(t)  =  t !/(*')  -  w0  ,  (8) 
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for  each  polygon  type. 


The  distortion  parameter  for  each  atom  site  complements  the  free  volume 
information.  The  pair  of  factors  of  free  volume  s/(i)  and  site  distortion 
parameter  u/(i)  are  indirect  structural  parallels  to  the  atomic  level  negative 
pressure  and  maximum  shear  stress. 

2.S  The  Structural  Relaxation  Simulation 


The  initial  state  for  the  present  simulation  of  structural  fluctuations  and 
relaxations  was  taken  to  be  the  state  of  the  material  in  the  melt  at  the  melting 
point  T *  =  0.25  for  the  2-D,  two-component  Cu  —  Zr  alloy  discussed  in  detail 
in  (I).  The  process  history  for  the  present  simulation  of  structural  alterations  is 
shown  in  Fig.  1,  in  which  two  complementary  structural  relaxations  have  been 
investigated:  (a)  a  relaxation  in  the  subcooled  liquid  region  between  Tm  and  Tg, 
at  a  temperature  of  T *  =  0.2,  and  (b)  a  relaxation  below  Tt  at  T*  =  0.1.  In 
both  of  these  cases,  the  material  is  suddenly  cooled  from  the  molten  state  at 
T*  =  0.25  (point  E  in  Fig.  1)  to  the  lower  temperature.  Thus,  the  cooling  rate 
to  the  isothermal  relaxation  temperatures  was  infinite  and  was  accomplished  by 
suddenly  decreasing  the  average  kinetic  energies  of  atoms  in  the  simulation  cell 
for  the  given  initial  coordinates  in  the  melt. 

We  discuss,  however,  also  in  considerable  detail,  the  nature  of  the  disorder 
and  its  continued  variation  at  T *  =  0.3  in  the  high-temperature  melt  to  obtain 
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a  better  appreciation  of  the  role  of  the  liquid-like  material  in  the  behavior  of  the 
melt.  All  simulations  were  carried  out  at  an  external  pressure  of  p*  =  1.0. 

m.  RESULTS  OF  THE  SIMULATION 

S.l  Structural  Defects  in  the  Amorphous  State 


In  (1),  we  have  demonstrated  that  in  the  amorphous  state  of  the  2-D  mate¬ 
rial  above  the  melting  point,  in  rapidly  solidified  form  as  a  sub-cooled  liquid, 
or  in  the  glassy  form  below  Tg,  the  disorder  that  is  evident  from  the  radial  dis¬ 
tribution  functions  (RDF)  of  atom  positions  is  not  uniformly  spread  over  the 
entire  material.  Instead,  there  is  a  considerable  fraction  of  the  material  which 
is  in  a  reasonably  well  ordered  form  at  any  snapshot  in  time,  being  made  up  of 
more  or  less  distorted  hexagonal  polygons.  In  the  melt,  these  ordered  regions 
were  found  separated  by  strings  of  “liquid-like”  material  percolating  through 
the  simulation  cell.  The  liquid-like  material  in  turn,  making  up  a  volume  frac¬ 
tion  of  about  0.4,  was  found  to  be  composed  of  primarily  polygons  of  5  and  7 
sides,  with  occasional  polygons  with  4  and  also  8  sides.  In  addition  to  this  de¬ 
fective  material  separating  more  ordered  domains,  occasional  dipoles  of  5  and  7 
sided  polygons  were  also  found  inside  the  quasi-ordered  domains.  Such  isolated 
dipoles  of  5-7  sided  polygons  in  the  quasi-ordered  domains  were  recognized  as 
edge  dislocations,  which  when  somewhat  displaced,  could  transform  into  two 
adjacent  pentagons  separated  by  a  flat  face  from  two  adjacent  large  distorted 
hexagons,  all  sharing  a  four  sided  corner.  This  configuration  was  identified  as 


11 


an  edge  dislocation  in  a  metastable  intermediate  condition  between  two  stable 
configurations  made  of  5-7  sided  polygons.  An  associated  study  of  the  mutual  ar¬ 
rangement  of  5-7  sided  polygon  pairs  established  that  most  of  such  arrangements 
where  two  5-7  sided  polygon  pairs  face  each  other  with  the  7  sided  polygons  be¬ 
ing  near  or  next  near  neighbors  and  the  5  sided  polygons  being  directed  outward, 
are  vacancies  or  vacancy  clusters.  On  the  other  hand,  the  arrangement  of  the  5-7 
sided  polygon  pairs  in  the  liquid-like  material  strings  separating  quasi-ordered 
domains  were  found  to  have  a  high  frequency  of  being  ordered  in  sequences  of 
touching.. ..5-7-5-7.... arrangements,  with  occasional  interruptions  breaking  this 
sequence.  Clearly,  in  a  perfect  hexagonal  mat,  such  regions  would  constitute 
high  angle  grain  boundaries,  and  in  the  amorphous  material,  they  must  be  in¬ 
terpreted  as  the  geometrically  necessary  arrangement  of  topological  defects  that 
result  in  the  misorientation  of  the  quasi-ordered  regions  which  they  separate. 
The  common  feature  of  all  of  these  defects  was  recognized  to  be  the  local  excess 
free  volume  which  they  provide,  which  was  established  to  be  on  the  average  a 
fraction  of  0.07  for  the  dipole  of  a  5-7  sided  polygon  pair. 

In  the  amorphous  material  of  the  subcooled  melt  above  Tt,  the  percolation 
condition  of  the  liquid-like  material  was  often  found  broken  in  snapshots  in 
time,  but  re-established  at  times  in  between  -  most  likely  due  to  the  decreasing 
misorientation  between  the  quasi-ordered  domains.  Below  the  glass  transition 
temperature,  the  percolation  condition  of  the  liquid-like  material  was  always 
found  broken. 
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S.2  Time  Dependent  Structural  Variations  in  the  Melt 


The  processes  of  structural  fluctuation  occurring  above  the  melting  point  at 
T *  =  0.3  were  followed  in  some  detail  over  a  period  of  500  time  steps  (25  thermal 
fluctuation  periods).  Figures  2a  -  2c  show  three  quite  representative  snapshots 
in  time,  after  300,  400,  and  500  time  steps.  It  is  clear  that  the  main  topological 
feature  is  the  more  or  less  aligned  strings  of  liquid-like  material  made  by  the 
association  of  5-7  sided  polygons,  with  a  dominant  alternating  sequence 
of.. ..5-7-5-7.... arrangements.  In  addition  to  the  presence  of  main  chains  of  strings 
of  polygon  pairs,  there  were  often  side  branches  either  joining  on  the  main  chain 
or  at  times  forming  isolated  subsidiary  chains.  In  addition,  however,  the  some¬ 
what  more  regular  hexagonal  regions  which  the  strings  separate  were  themselves 
quite  markedly  distorted.  As  discussed  in  (I),  at  this  state  of  the  material,  the 
2-D  RDF  shows  only  3  broad  and  interconnected  peaks,  made  up  of  the  fusion 
of  the  second  and  third  peaks,  and  the  fourth  and  fifth  peaks,  in  addition  to 
the  first  peak.  A  certain  number  of  4  edged  corners  separating  two  neighbor¬ 
ing  5  sided  polygons  were  visible  in  all  the  frames.  As  discussed  above,  these 
are  higher  energy  configurations  of  the  elementary  edge  dislocation  cores  in  the 
idealization  of  the  liquid-like  regions  as  distorted  high  angle  grain  boundaries. 
Since  the  initial  state  of  the  2-D  material  before  being  melted  consisted  of  a  field 
of  perfect  hexagons  subjected  to  periodic  boundary  conditions,  certain  conserva¬ 
tion  laws  must  apply  to  the  defects  in  the  disordered  states  of  the  material  in  the 
melt  or  below  the  melting  temperature.  These  conservation  laws  were  studied 
in  considerable  detail  by  Smith  [17]  in  2-D  area  filling  polygons.  In  one  limiting 
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form,  where  only  5,  6,  and  7  sided  polygons  exist  in  the  field,  the  numbers  Pn 
of  the  n  sided  polygons  relate  to  the  numbers  C  of  the  4  edged  corners  by  the 
conservation  law  of: 


£(6  -  n)Pn  -  2C  =  0  .  (9) 


Comparison  of  the  three  figures  of  Fig.  2a-2c  shows  that  in  the  melt,  the 
liquid-like  material  of  the  5-7  sided  polygon  regions  do  not  remain  stationary  in 
time,  but  are  constantly  shuffled,  maintaining,  however,  a  percolation  condition 
on  the  average.  Even  so,  percolating  chains  often  break  up  and  establish  asso¬ 
ciations  with  other  subsidiary  chains,  as  apparently  the  quasi-ordered  domains 
shift  and  rotate  randomly.  In  this  manner,  the  liquid-like  material  is  capable 
of  wandering  over  the  cubic  2-D  field,  and  can  in  principle  have  access  to  every 
point  in  the  2-D  space.  It  is  useful  to  consider  this  fluctuating  disorder  as  a 
graphical  representation  of  a  configurational  entropy  se  attributable  to  the  pres¬ 
ence  of  the  liquid-like  material  alone,  which  could  be  defined  along  conventional 
lines  as: 


se  =  ktn 


(  W  2)1  \ 

V((AT/2) -n)!n!j  ’ 


(10) 
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where  N  is  the  total  number  of  atoms  in  the  cell,  and  n  is  the  number  of  the  5-7 
sided  pairs  of  polygons.  The  maximum  in  this  configurational  entropy  occurs 
when  n  =  N/4,  or  at  a  value  of  the  fraction  of  the  liquid-like  material  of  0.5. 
The  actual  fractions  of  liquid-like  material  made  up  of  the  5-7  sided  polygon 
pairs  at  the  melting  point  was  found  to  vary  between  a  value  of  0.36  to  a  high 
value  of  0.41.  On  the  basis  of  this,  the  configurational  entropy  is  found  to  have 
reached  95%  of  its  peak  value  given  by  the  expression  in  Eqn.  (10).  Clearly,  in 
the  free  energy  of  the  liquid  state,  the  internal  potential  energy  of  the  atoms  also 
plays  an  important  role,  and  limits  the  role  of  the  configurational  entropy.  This 
is  manifested  by  the  quasi-ordered  regions  which  force  the  disordered  liquid-like 
material  into  the  border  regions  between  the  ordered  regions.  Nevertheless,  the 
expression  in  Eqn.  (10)  furnishes  a  good  guide  to  account  for  the  overall  fraction 
of  the  liquid-like  material  in  the  melt. 

Examination  of  the  break-up  and  re-formation  of  the  percolating  strings  of 
the  5  and  7  sided  polygons  has  shown  that  when  such  breaks  occur  through 
the  transformation  of  these  polygons  into  distorted  hexagons,  usually  the  "dis¬ 
solved  misfit"  re-appears  in  the  immediate  neighborhood  as  a  pair  of  new  5-7 
sided  polygons,  associating  themselves  to  the  ends  of  other  neighboring  chains. 
This  is  to  be  expected  in  the  crystal  defect  idealization,  where  terminations  of 
chains  of  5-7  sided  polygons  would  represent  terminations  of  high  angle  grain 
boundaries  with  a  relatively  high  stress  that  would  attract  other  5-7  sided  poly¬ 
gons  to  provide  compensation  of  misfit.  Alternatively,  viewing  the  boundaries 
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of  5-7  sided  polygons  as  high  angle  grain  boundaries  separating  quasi-ordered 
domains,  shifts  in  the  domains  will  move  the  boundaries  around,  including  direc¬ 
tions  where  the  geometrically  necessary  arrangement  of  dislocation  cores  requires 
occasional  gaps  between  the  5-7  sided  polygon  dipoles. 

In  a  more  elementary  form,  gradual  drifts  of  isolated  5-7  sided  polygon 
dipoles  towards  each  other,  have  also  been  observed  occasionally.  Such  events, 
however,  are  quite  rare,  since  isolated  5-7  sided  polygon  dipoles  are  themselves 
rare,  at  or  around  the  melting  point,  where  their  combined  glide  and  climb 
mobility  could  make  a  contribution  to  structural  fluctuations. 

Parenthetically,  at  low  temperatures  in  well-relaxed  structures  of  a  high 
crystalline  content,  where  such  isolated  5-7  polygon  pairs  occasionally  occur  and 
stand  out,  their  mobility  is  quite  low.  Nevertheless,  in  the  structural  relaxation 
process  leading  to  material  with  increasing  perfection  through  the  disappear¬ 
ance  of  the  liquid-like  material,  such  pair-wise  eliminations  of  5-7  sided  dipoles 
is  topologically  necessary,  as  we  will  demonstrate  below.  One  such  case  of  glide- 
controlled  approach  of  two  edge  dislocations  toward  each  other  under  their  mu¬ 
tual  attraction  in  a  partially  relaxed  glass  at  T*  =0.1  in  the  course  of  further 
relaxation  is  shown  in  the  sequence  of  Fig.  3a-3e,  in  the  upper  right  corner  of 
the  figures  marked  with  arrows.  The  figure  also  shows  clearly  two  metastable 
configurations  in  the  slowly  gliding  edge  dislocation  as  two  adjoining  polygons 
with  5  edges  sharing  a  4  edged  corner.  The  dislocations  are  marked  in  Fig.  3a 
with  the  appropriate  Vs,  to  identify  them  and  their  mutual  arrangement  that 
should  result  in  an  attraction.  Parenthetically,  examinations  of  specific  configu¬ 
rations  in  well-relaxed  material  at  T*  =  0.1  have  shown  several  instances  where 
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an  isolated  edge  dislocation  has  assumed  the  metastable  configuration  of  two  5 
sided  adjacent  polygons  only  to  return  to  its  initial  configuration  of  a  pair  of  5-7 
sided  polygons.  This  indicates  that  the  two  neighboring  5-7  sided  polygon 
configurations  are  the  lower  energy  form  and  that  the  two  5  sided  polygons 
indeed  represent  a  metastable  form  of  somewhat  higher  energy. 

S.8  Structural  Relaxation  in  the  Subcooled  Melt,  Tt  <  T  <  Tm 


The  sequence  of  structural  relaxations  at  a  temperature  of  7"  =  0.2  between 
Tt  and  Tm  in  the  subcooled  melt  is  shown  in  Figs.  4a-4d,  starting  at  a  time  of 
27  fluctuation  periods,  after  initiation  of  the  isothermal  relaxation,  Fig.  4a,  and 
continuing  to  100  fluctuation  periods,  Fig.  4d.  Over  this  period  of  structural 
relaxation,  it  is  evident  that  the  overall  content  of  the  liquid-like  material  is  de¬ 
creasing  monotonically,  as  it  is  expected  to  do.  The  overall  fraction  of  liquid-like 
material  in  Fig.  4a,  after  27  fluctuation  periods  has  decreased  from  the  starting 
fraction  of  about  0.4  in  the  melt  to  about  0.31.  In  this  range,  the  percolating 
strings  of  liquid-like  material  are  almost  intact,  up  to  a  total  relaxation  period 
of  35.  This  is  roughly  the  level  of  structural  alterations  studied  in  (I),  where 
Tt  was  established  at  a  cooling  rate  involving  a  time  period  somewhat  less  than 
this  relaxation  time.  In  Fig.  4d,  after  a  total  relaxation  time  of  100  fluctuation 
periods,  the  total  fraction  of  liquid-like  material  has  decreased  down  to  0.19,  and 
the  percolation  condition  has  clearly  been  broken,  as  the  mutual  misorientation 
between  quasi-ordered  domains  has  very  markedly  decreased.  Thus,  this  end 
state  of  the  material,  associated  with  a  slower  rate  of  cooling,  is  now  no  longer 
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a  subcooled  liquid,  but  a  well-relaxed  glass.  Apart  from  the  gradual  reduction 
of  the  liquid-like  fraction  of  the  material,  the  kinematics  of  the  structural  relax¬ 
ation  at  T *  =  0.2,  in  the  range  of  subcooled  liquid,  exhibit  most  of  the  modes 
of  structural  variations  found  in  the  high  temperature  melt. 

8.4  Structural  Relaxation  below  Tt,  at  T*  =  0.1 


Figures  5a-5e  give  a  sequence  of  changes  in  the  distribution  of  liquid-like  ma¬ 
terial  over  a  time  of  200  fluctuation  periods  at  T*  =0.1.  While  at  the  beginning, 
in  the  as-quenched  structure,  a  condition  of  near  percolation  exists,  this  disap¬ 
pears  quite  rapidly.  Over  the  period  of  simulation,  the  fraction  of  the  liquid-like 
material  decreases  monotonically  from  0.33  to  0.049.  Since  the  liquid-like  mate¬ 
rial  is  made  up  of  structural  dipoles  of  5-7  sided  polygons,  which  should  usually 
attract  each  other  and  result  in  cancellation  when  such  pairs  in  exact  compen¬ 
sating  configurations  encounter  each  other.  This  cancellation  might  be  expected 
to  obey  second  order  kinetics,  i.e., 

t  =  •  <“> 


where  /  is  the  fraction  of  liquid-like  material,  and  a  is  primarily  the  reciprocal 
of  a  characteristic  time  constant.  If  this  were  indeed  the  case,  a  specific  form  of 
time  dependent  decrease  should  result,  that  should  be: 
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l/f-l/fo  =  at  . 


(12) 


That  this  is  indeed  the  case,  is  shown  in  Fig.  6.  From  the  figure,  we  establish 
that  the  characteristic  time  constant  1/a  =  15  fluctuation  periods.  This  gives  a 
first-order  fit  to  the  actually  more  complex  and  distributed  kinetics  of  structural 
relaxation  that  will  be  the  subject  of  a  separate  communication  [13]  (referred 
to  earlier  as  (III)).  Ahn  and  Li  [18 j  have  interpreted  the  kinetics  of  the  strain 
recovery  experiments  of  Berry  [19]  in  PdgoSt20,  with  second  order  kinetics  on 
the  basis  of  an  analog  to  cancellation  of  dislocation  pairs.  Our  simulation  is  in 
agreement  with  their  findings. 

The  steps  in  the  elimination  of  liquid-like  material  are  usually  complex. 
Nevertheless,  in  the  final  stages  of  this  elimination  where  a  limiting  idealization 
of  such  defects  as  edge  dislocation  cores  in  boundaries  is  useful,  this  point  of 
view  clarifies  the  kinetics  of  these  processes  considerably.  Two  examples  of  mu¬ 
tual  cancellation  of  5-7  sided  polygon  pairs  that  have  been  observed  are  shown 
in  Fig.  7a  and  7b.  The  pair  in  Fig.  7a  acts  like  two  slightly  displaced  opposite 
edge  dislocations,  which  are  not  exactly  on  the  same  plane  to  cancel  fully.  The 
accompanying  symbolic  sketch  below  the  polygons  suggests  that  the  two  pairs 
attract  to  cancel.  In  the  intermediate  figure,  they  do  not  produce  full  cancella¬ 
tion,  but  result  in  a  momentary  partial  interstitial,  which  in  the  last  frame,  is 
dispersed  among  the  four  hexagons  by  a  small  rearrangement  of  polygons.  In 


19 


Fig.  7b,  the  two  approaching  dislocations  have  Burgers  vectors  making  a  60° 
angle  with  each  other,  and  are  in  a  mutually  attracting  configuration. 

Upon  contact,  they  react  to  give  a  single  product  dislocation  with  a  Burgers  vec¬ 
tor  making  a  60°  angle  with  both  of  the  previous  dislocations,  as  shown  in  the 
accompanying  symbolic  sketch.  Clearly,  the  sequences  of  the  relaxation  shown 
in  Figs.  5a-5e  take  the  structure  from  a  well-relaxed  glass  progressively  to  a 
crystal.  On  such  a  small  simulation  plane,  the  transition  is  gradual,  and  its  on¬ 
set  cannot  be  clearly  identified.  This  difficulty,  however,  is  common  in  the  study 
of  atomic  amorphous  media,  where  no  single  threshold  criterion  has  proved  fully 
satisfactory. 

The  results  of  the  structural  relaxations  at  T*  =0.1  over  200  fluctuation 
periods  (4000  time  steps)  are  typical  of  all  structural  relaxations,  which  we 
have  studied  and  are  most  instructively  viewed  by  comparing  the  changes  in  the 
distribution  of  a  number  of  intensive  local  atomic  properties,  before  and  after 
the  entire  relaxation  period.  These  results  for  the  average  atomic  level  enthalpy, 
Voronoi  polygon  volume  (an  indirect  measure  of  free  volume),  site  distortion 
parameter,  atomic  level  pressure,  atomic  level  maximum  shear  stress,  and  atomic 
level  bulk  modulus,  are  shown  in  Figs.  8-13,  first  as  averages  over  all  atoms  (a), 
then  for  averages  over  only  the  softer  Cu  atoms  (b),  and  finally,  for  averages  over 
only  the  stiffer  Zr  atoms  (c).  In  the  figures,  only  the  changes  in  the  distributions 
have  been  shown  that  modify  the  initial  distributions  into  the  final  distributions. 
In  the  cases  of  atomic  levels  of  enthalpy,  Voronoi  polygon  volume,  and  site 
distortion  parameter,  the  extremes  of  excess  properties  have  been  markedly 
diminished  to  enrich  the  distributions  at  the  low  end.  This  is  particularly  clear 
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for  the  enthalpies  and  site  distortion  parameters,  while  in  the  case  of  the  Voronoi 
polygon  volume,  both  the  large  and  the  small  VP  volumes  have  been  eliminated 
in  favor  of  the  intermediate  volume  VPs,  i.e.,  the  hexagons.  There  are  significant 
differences  between  the  Zr  atom  with  its  stiffer  environment,  and  the  Cu  atom 
with  its  softer  environment.  The  structural  misfit  produced  by  the  disorder  is 
more  readily  accommodated  by  the  softer  Cu  atoms,  as  would  be  expected  from 
the  analogous  case  of  a  random  network  of  interconnected  stiff  and  compliant 
springs.  In  such  a  network  with  incorporated  random  misfit,  more  energy  is 
stored  in  the  more  compliant  springs  than  in  the  stiffer  springs,  i.e.,  in  the  Cu 
sites,  rather  than  in  the  Zr  sites.  Thus,  upon  structural  relaxation,  the  changes 
in  excess  enthalpy  are  larger  for  the  Cu  sites.  Equally  dramatic  differences  are 
apparent  in  the  changing  atomic  level  pressure  distributions  between  the  Cu  and 
the  Zr  sites.  Evidently,  in  the  initial  highly  disordered  glassy  state,  the  stiffer 
and  larger  Zr  atoms  form  an  interconnected  super  structure  in  compression 
with  the  interspersed  Cu  atoms  being  generally  under  negative  pressure.  The 
additional  applied  pressure  in  the  simulation  is  also  largely  born  by  the  Zr  super 
structure.  Upon  structural  relaxation,  steadily  better  fits  among  the  atoms 
are  achieved,  no  doubt  because  of  an  increasing  chemical  short  range  order, 
as  has  been  already  reported  by  Kabayashi  and  Takeuchi  [10].  As  a  result, 
the  Cu  atoms  are  subjected  to  less  negative  pressure,  and  the  Zr  atoms  are 
subjected  to  less  positive  pressure.  Hence,  the  incremental  atomic  level  pressure 
is  positive  on  the  Cu  atoms,  and  negative  on  the  Zr  atoms.  The  decrease  of  the 
negative  pressure  from  the  Cu  atoms  is  also  evident  from  an  increase  in  their 
bulk  modulus. 
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rv.  DISCUSSION 


4.1  Characteristics  of  the  5  and  7  aided  Polygon  Pairs 


In  (I),  we  have  demonstrated  that  in  a  distorted  field  of  area  filling  polygons, 
the  average  heptagons  have  a  10%  larger  area  than  the  average  hexagons,  while 
the  average  pentagons  are  only  3%  smaller  than  the  latter  -  making  the  pair 
of  pentagon  and  heptagon,  i.e.,  the  5-7  pair,  have  an  excess  area  of  7%  over 
two  average  hexagons.  Considering  that  the  average  hexagon  at  the  melting 
point  itself  has  an  average  area  that  is  about  9%  larger  than  its  area  of  0°K , 
the  5-7  sided  pair  of  polygons  are  16%  larger  than  two  reference  hexagons  at 
rest.  Such  an  areal  expansion  should  produce  significant  loss  of  cohesion  in  any 
realistic  inter-atomic  binding  potential  and  therefore,  constitutes  an  element  of 
liquid-like  material.  In  the  polygon  mats  of  molten  material,  occasional  octagons 
and  squares  were  also  found.  These,  however,  were  always  part  of  the  string  of 
liquid-like  material,  and  were  interpreted  as  further  distortions  of  7  sided  and  5 
sided  polygons. 

Since  the  5-7  sided  polygon  pairs  were  found  to  hold  the  key  to  many  struc¬ 
tural  and  mechanical  relaxation  phenomena,  some  of  their  other  properties  were 
also  evaluated.  In  Fig.  14a,  the  change  in  the  site  distortion  parameter  is  shown 
with  the  changing  number  of  edges  of  polygons  in  a  simulation  at  T*  =0.1.  The 
solid  line  gives  the  change  in  the  raw  measure  of  the  site  distortion,  as  obtained 
from  an  evaluation  of  Eqn.  (7).  If  corrections  are  made  according  to  Eqn.  (8), 
however,  subtracting  the  trivial  decrease  in  the  distortion  parameter  with 
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increasing  numbers  of  sides  of  polygons  given  in  Table  I,  the  information  shown 
by  the  broken  lines  in  Fig.  14a  is  obtained,  which  indicates  that  the  average  pair 
of  hexagons  are  far  less  distorted  than  a  5-7  pair  of  polygons.  The  corresponding 
evaluation  for  the  average  excess  enthalpy  of  polygons  over  that  for  the  average 
hexagon,  as  evaluated  by  Eqn.  (6),  is  shown  in  Fig.  14b.  Thus  it  is  clear  that 
the  excess  free  volume  in  the  5-7  sided  polygon  pairs  results  in  excess  enthalpy. 

4.8  Topological  Features  of  Structural  Relaxation 


The  surprising  aspect  of  the  two-dimensional  simulation  of  structural  fluc¬ 
tuations  in  the  melt  and  the  features  of  its  relaxations  below  the  melting  point 
are  the  striking  parallels  that  it  has  revealed  to  a  defective  crystalline  material. 
In  our  simulation,  a  melt  is  made  up  of  very  small  quasi-ordered  domains  of 
more  or  less  distorted  hexagons  separated  by  boundaries  of  much  more  severely 
distorted  polygons,  which  are  analogous  to  high  angle  grain  boundaries.  This 
distorted  boundary  material  was  found  to  be  predominantly  made  up  of  pairs 
of  polygons  with  5  and  7  sides  arranged  mostly  in....5-7-5-7....  sequences,  with 
occasional  breaks  bridged  by  a  severely  distorted  hexagon.  In  the  melt,  this 
boundary  material  was  found  to  occupy  about  40%  of  the  total  area.  Since  the 
boundary  material  is  on  the  average  one  atom  thick,  the  quasi-ordered  domains 
that  it  separates  are  only  3-5  atoms  across.  Thus,  the  quasi-ordered  domains 
should  not  be  viewed  as  individual  grains,  but  at  best,  only  as  topologically 
short  range  ordered  material.  Since  the  polygon  pairs  with  5  and  7  sides  rep¬ 
resent  edge  dislocations  in  a  perfect  mat,  the  boundary  material  where  they 
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are  sequentially  arranged  must  be  viewed  in  abstraction  as  a  high  angle  grain 
boundary  made  up  entirely  of  dislocation  cores,  having  collectively  no  long  range 
stress  field  characteristic  of  dislocations.  In  fact,  the  only  dominant  property  of 
the  boundary  material  with  the  5  and  7  sided  polygons  is  an  excess  volume  (free 
volume)  of  about  7%  larger  than  the  average  volume  of  the  distorted  hexagons  in 
the  quasi-ordered  domains,  which  themselves  are  about  9 %  larger  at  the  melting 
point  than  the  reference  hexagons  of  a  perfect  mat  at  Q°K.  Thus,  the  boundary 
material  has  an  average  dilatation  of  about  16%  in  comparison  with  a  perfect 
mat  at  0°K.  This  should  reduce  the  inter-atomic  cohesion  in  the  boundaries 
to  a  relatively  low  level.  Following  Cohen  and  Grest  [20),  we  have  called  the 
boundary  material  “liquid-like” ,  and  have  noted  that  it  percolates  through  the 
entire  2-D  space.  That  this  material  is  indeed  liquid-like,  and  that  it  prescribes 
the  sites  for  shear  transformations  when  the  mats  are  sheared,  is  presented  and 
discussed  in  great  detail  in  the  last  of  the  accompanying  communications,  which 
we  will  refer  to  as  IV  [21]. 

Occasional  edge  dislocations  are  also  found  in  quasi-ordered  domains  which 
have  shown  mutual  attractions  to  each  other,  have  moved  by  glide,  and  in  a  few 
instances,  fused  with  other  dislocations  to  give  cancellation  of  pairs  or  expected 
reaction  products.  Such  events,  however,  where  isolated  well-delineated 
dislocations  have  given  purely  dislocation  related  glide,  and  associated  phenom¬ 
ena  were  quite  rare.  Nevertheless,  in  the  structural  relaxations,  increasing  order 
and  decreasing  misorientation  in  the  quasi-ordered  domains  through  the  elim¬ 
ination  of  the  5-7  sided  pairs  of  polygons  as  dislocation  cores  in  the  boundary 
material,  could  only  have  happened  by  pair-wise  cancellations  of  the  5-7  sided 
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polygon  pairs.  The  reasonably  good  fit  to  second  order  kinetics  of  the  reduction 
of  the  fraction  of  boundary  material  is  in  support  of  this  sequence  of  events. 


The  results  of  our  simulation  are  in  remarkable  agreement  with  the  observa¬ 
tions  of  Fukushima  and  Ookawa  [22],  who  carried  out  quite  detailed  pioneering 
experiments  on  melting  with  vibrating  Bragg  soap  bubble  rafts.  They  reported 
that  in  their  rafts,  melting  occurred  when  a  high  concentration  of  high  angle 
grain  boundaries  divided  the  raft  up  into  small  isolated  quasi-crystalline 
domains  of  roughly  the  same  dimensions  as  our  quasi-crystalline  domains.  Since 
the  inter-bubble  potential  provides  a  good  match  to  interatomic  pair  potentials 
for  close-packed  metal  atoms  in  the  binding  region  (for  a  comparison  of  the 
potentials,  see  Shi  and  Argon  [23]}  their  findings  should  not  be  surprising  - 
even  though  the  vibrational  amplitudes  of  bubbles  must  have  had  a  substantial 
component  out  of  the  plane  of  the  raft. 

4.S  Generalization  of  Results  to  Three  Dimensions 


Clearly,  the  results  of  our  simulation  must  be  taken  with  some  caution.  As 
already  mentioned  in  (I),  two-dimensional  disordered  media  entraps  larger  free 
volume  than  three-dimensional  matter  -  if  the  disorder  were  indeed  uniformly 
spread  out.  It  can  be  argued,  however,  that  the  partitioning  of  order  and  dis¬ 
order  could  not  be  much  different  in  three-dimensional  amorphous  media.  Cer¬ 
tainly,  the  make  up  of  the  percolating  liquid-like  boundary  surfaces  in  three- 
dimensional  material  will  have  a  larger  collection  of  polyhedra  that  provides 
a  transition  between  the  different  quasi-lattices  of  the  small  ordered  domains. 


25 


Self-consistent  model  studies  of  packing  of  rigid  spheres  randomly  in  3-D 
spaces  have  indicated  that  the  mechanical  properties  of  the  medium  assume 
rigid  behavior  when  a  percolation  criterion  is  reached  at  a  volume  fraction  of 
0.4  [24,25],  which  is  almost  identical  to  the  value  for  the  liquid-like  material 
boundaries  in  our  complementary  case  of  2-D  simulation,  where  this  fraction  of 
material  is  deformable.  Although  the  spatial  arrangement  of  the  quasi-ordered 
regions  in  3-D  space  cannot  be  determined  from  these  studies,  it  is  quite  likely 
that  they  are  equiaxed  domains  leaving  behind  walls  of  a  closed-cell  network 
or  liquid-like  material,  rather  than  being  arranged  in  a  spoke  network.  Thus, 
we  expect  that  most  of  the  qualitative  observations  from  our  2-D  simulation 
should  be  applicable  to  3-D  materials,  particularly  on  topological  features  and 
on  the  inelastic  mechanical  response.  The  latter  will  be  the  subject  of  our  last 
communication  (IV)  [21]. 
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TABLE  I 


Base  Levels  of  Distortion  Parameter  w0  for  regular  Polygons 


Polygon 

1V0 

Equilateral  triangle 

1.286 

Square 

1.128 

Pentagon 

fllla 

Hexagon 

a 

Heptagon 

Octagon 

BEa 
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FIGURE  CAPTIONS 


Fig.  1  Schematic  representation  of  quenching  history  in  the  structural  relax¬ 
ation  simulations:  ABCD  heating  path  leading  to  melting  and  super¬ 
heating,  DE  slow  cooling  in  equilibrium,  EF  cooling  at  infinite  rate, 
and  FG  isothermal  structural  relaxation. 

Fig.  2  Structural  fluctuations  in  a  high  temperature  melt  at 

T *  =  0.3  (=  1.2  Tm):  (a)  after  15  fluctuation  periods  rj  (b)  after  20 
fluctuations,  (c)  after  25  fluctuations.  Note  the  distribution  of  liquid¬ 
like  material  made  up  of  5-7  sided  polygon  pairs. 

Fig.  3  A  sequence  of  glide-like  approach  of  two  edge  dislocations  in  a  quasi- 
ordered  domain  at  T*  =  0.1.  The  position  and  nature  of  dislocations 
are  indicated  with  arrows  and  letters  T  in  frame  (a).  Frames  b-e  show 
continued  stages  of  approach. 

Fig.  4  Structural  relaxation  sequences  in  a  subcooled  melt  at  T *  =  0.2 
(Tt  <  T  <  Tm):  (a)  after  t  =  27 *7,  (b)  t  =  30*7, 

(c)  t  =  35 *7,  and  (d)  t  =  100*7,  showing  progressive  loss  of  liquid-like 
material. 

Fig.  5  Structural  relaxation  sequences  in  a  glass  at  T *  =  0.1  (T  <  Tt ): 

(a)  t  —  0,  (b)  t  =  20tj,  (c)  t  =  70>;,  (d)<  =  120*7,  and  (e)  t  =  200*7. 

Fig.  6  Time  dependent  reduction  in  fraction  of  liquid-like  material  at 

T *  =  0.1  during  monotonic  structural  relaxation.  The  relationship 
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indicates  second  order  kinetics  arising  from  pair-wise  cancellation  of 
polygon  dipoles. 

Fig.  7  Two  sequences  of  cancellation  of  dipoles:  (a)  dipole  cancellation 

analogous  to  two  opposite  edge  dislocation  cores,  (b)  dipole  reduction 
analogous  to  reaction  of  two  60°  edge  dislocation  cores. 

Fig.  8  Total  incremental  effect  of  structural  relaxation  at  T*  =  0.1  on  atomic 
level  enthalpy:  (a)  all  atoms,  (b)  Cu  atoms  only,  and  (c)  Zr  atoms 
only. 

Fig.  9  Total  incremental  effect  of  structural  relaxation  at  T*  —  0.1  on  volume 
per  atom  (Voronoi  polygon  volume):  (a)  all  atoms,  (b)  Cu  atoms 
only,  and  (c)  Zr  atoms  only. 

Fig.  10  Total  incremental  effect  of  structural  relaxation  at  T*  —  0.1  on  atomic 
site  distortion  parameter:  (a)  all  atoms,  (b)  Cu  atoms  only, 
and  (c)  Zr  atoms  only. 

Fig.  11  Total  incremental  effect  of  structural  relaxation  at  T*  =  0.1  on  atomic 
level  pressure:  (a)  all  atoms,  (b)  Cu  atoms  only,  and  (c)  Zr  atoms 
only. 

Fig.  12  Total  incremental  effect  of  structural  relaxation  at  T*  =  0.1  on  atomic 
level  maximum  shear  stress:  (a)  all  atoms,  (b)  Cu  atoms  only,  and 
(c)  Zr  atoms  only. 

Fig.  13  Total  incremental  effect  of  structural  relaxation  at  T'  =  0.1  on  atomic 
level  bulk  modulus:  (a)  all  atoms,  (b)  Cu  atoms  only,  and 
(c)  Zr  atoms  only. 
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Fig.  14  Comparison  of  intensive  properties  of  polygons  with  six,  five,  and 
seven  sides  for  partially  relaxed  mats  at  T*  =0.1:  (a)  corrected 
distortion  parameter  for  average  hexagon  is  less  than  those  for 
average  pentagon  and  heptagon,  (b)  atomic  level  enthalpy  is  lower 
for  hexagons  than  for  pentagons  and  heptagons. 
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-  Schematic  representation  of  quenching  history  in  the  structural 
relaxation  simulations:  ABCD heating  path  leading  to  melting  and 
superheating,  DE  slow  cooling  in  equilibrium,  EF  cooling  at  in¬ 
finite  rate,  and  FG  isothermal  structural  relaxation. 


Structural  fluctuations  in  a  superheated  melt  at  T  =*  0.3 
(*  1*2  Tm):  (a)  after  15  fluctuation  periods  t,  (b)  after  20 
fluctuations,  (c)  after  25  fluctuations.  Note  the  distribution 
of  liquid-like  material  made  up  of  5-7  sided  polygon  pairs. 


Fig.  2b 


A  sequence  of  glide-like  approach  of  two  edge  dislocations  in  a 
quasi-ordered  domain  at  T*s  0.1.  The  position  and  nature  of 
dislocations  are  indicated  with  arrows  and  letters  T  in  frame 
(a).  Frames  b-e  show  continued  stages  of  approach. 


Fig.  3b 


Fig.  3d 


-  Structural  relaxation  sequences  in  a  subcooled  melt  at  T  «  0.2 
(Tg  <  T  <  Tm):  (a)  after  t  *  27t,  (b)  t  ■  30x,  (c)  t  *  35t,  and 
(dj  t  =  IOQt,  showing  progressive  loss  of  liquid-like  material. 
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Fig.  5  -  Structural  relaxation  sequences  in  a  glass  at  T  *  0.1  (T  <  Tj 
(a)  t  *  0,  (b)  t  -  20t,  (c)  t  *  70t,  (d)  t  ■  120t,  and  9 
(e)  t  3  200t. 


Fig.  5b 


Fig.  5c 


Fig.  6  -  Time  dependent  reduction  in  fraction  of  liquid-like  material  at 
T  =  0.1  during  monotonic  structural  relaxation.  The  relation¬ 
ship  indicates  second  order  kinetics  arising  from  pair-wise 
cancel letion  of  polygon  dipoles. 
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Fig.  7 


Two  sequences  of  cancellation  of  dipoles:  (a)  dipole  cancella- 
tion  analogous  to  two  opposite  edge  dislocation  cores, 

(b;  dipole  reduction  analogous  to  reaction  of  two  60°  edqe  dis¬ 
location  cores.  3 
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-  Total  incremental  effect  of  structural  relaxation  at  T  =  0.1 
atomic  level  enthalpy:  (a)  all  atoms,  (b)  Cu  atoms  only,  and 
(c)  Zr  atoms  only. 
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-  Total  incremental  effect  of  structural  relaxation  at  T  =  0.1  on 
volume  per  atom  (Voronoi  polygon  volume):  (a)  all  atoms, 

(b)  Cu  atoms  only,  and  (c)  Zr  atoms  only. 
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Fig.  10  -  Total  incremental  effect  of  structural  relaxation  at  T 
atomic  site  distortion  parameter:  (a)  all  atoms,  (b)  Cu 
only,  and  (c)  Zr  atoms  only. 
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Total  incremental  effect  of  structural  relaxation  at  T*=  0.1 
atomic  level  pressure:  (a)  all  atoms,  (b)  Cu  atoms  only,  and 
(c)  Zr  atoms  only. 
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LOCAL  SHEAR  STRESS  (A  ATOMS) 
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£lSL»  ^2  -  Total  incremental  effect  of  structural  relaxation  at  T  *  0.1  on 
atomic  level  maximum  shear  stress:  (a)  all  atoms,  (b)  Cu  atoms 
only,  and  (c)  Zr  atoms  only. 
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LOCAL  BULK  MODULUS  (B  ATOMS) 
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Total  incremental  effect  of  structural  relaxation  at  T*=  0.1  i 
atomic  level  bulk  modulus:  (a)  all  atoms,  (b)  Cu  atoms  only, 
and  (c)  Zr  atoms  only. 
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-  Comparison  of  intensive  properties  of  polygons  with  six,  five, 
and  seven  sides  for  partially  relaxed  mats  at  T  *  0.1: 

(a)  corrected  distortion  parameter  for  average  hexagon  is  less 
than  those  for  average  pentagon  and  heptagon,  (b)  atomic  level 
enthalpy  is  lower  for  hexagons  than  for  pentagons  and  heptagons. 
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